Can Planck-scale physics be seen in the cosmic microwave background ? 
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We investigate the potential of observations of anisotropies in the cosmic microwave background 
(CMB) and large-scale structure in the Universe to detect possible modifications of standard in- 
flationary models by physics beyond the Planck scale. A generic model of the primordial density 
fluctuations is investigated, and we derive constraints on its parameters from current data. We 
conclude that the currently available data do not put very stringent constraints on this model. Fur- 
thermore, we use simulated power spectra from the Sloan Digital Sky Survey (SDSS) and the Planck 
satellite to show that it is unlikely that a trans-Planckian signature of this type can be detected in 
CMB and large-scale structure data. 

PACS numbers: 04.60.-m, 98.80.-k, 98.80.Cq 
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I. INTRODUCTION 



Quantum gravitational effects are expected to mod- 
ify the spectrum of primordial density fluctuations pro- 
duced during the inflationary phase in the very early Uni- 
verse. Since the primordial fluctuations are the seeds 
for the anisotropies in the cosmic microwave background 
(CMB) and for the large-scale structures we observe in 
the Universe, cosmological observations have the poten- 
tial to shed light on on Planck-scale physics. However, 
although enormous progress has been made in observa- 
tional cosmology in the past few years, it is a highly 
non-trivial task to separate the primordial density fluc- 
tuations from the present-day, processed power spectra 
which we can measure. In this paper we show that 
the currently available data do not impose any strin- 
gent constraints on trans-Planckian modulations of the 
primordial power spectrum. Furthermore, we investigate 
whether future, more precise measurements will improve 
the prospects for detecting a signature for quantum grav- 
ity. Based on suggestions for the form and size of trans- 
Planckian effects on the primordial power spectrum (see 
e-g- Q> IE IE ^ IE Hi Q)' consider a model where a 
power-law spectrum is modulated byan oscillating func- 
tion of the comoving wavenumber [EISH- perform- 
ing a full likelihood analysis of simulated data sets, we 
are led to the conclusion that it will be extremely difficult 
to detect trans-Planckian effects in CMB and large-scale 
structure data. 

The structure of this paper is as follows. In section II 
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we give a brief summary of the work on trans-Planckian 
effects on the primordial power spectrum and introduce 
our model. Section III contains a discussion of the cos- 
mological data, with an emphasis on the unavoidable 
smearing of power by the window functions. Section IV 
describes our likelihood analysis and the results. 



II. TRANS-PLANCKIAN PHYSICS AND THE 
PRIMORDIAL POWER SPECTRUM 

The standard calculation of the perturbations pro- 
duced during the inflationary phase is based on flat-space 
quantum field theory, and the initial conditions for the 
inflaton (j) are imposed in the infinite past. Since the 
Universe is expanding, this means that any given Fourier 
mode k was infinitesimally small at that time. Since our 
expectation is that 'new physics' becomes important on 
length scales below the Planck length, the standard cal- 
culation is probably too naive, and we should expect 
corrections from quantum gravitational effects. We do 
not know what the correct fundamental theory on scales 
comparable to and smaller than the Planck length is, 
and hence it is impossible to predict these corrections 
with a high degree of certainty. One can only devise 
plausible scenarios based on informed guesses as to what 
the significant features of a theory of quantum gravity 
should be like. Many proposals have been made with 
a varying degree of optimism regarding possible observ- 
able signatures 0, ji, J,i ^ IE ■ important point 
is whether the leading-order correction to the primordial 
power spectrum is linear in H/A, where H is the Hub- 
ble constant during inflation and A is the energy scale 
where new physics enters. If the correction is of higher 
order than linear, it will certainly be undetectable, so we 
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will assume that it is linear. With one specific choice of 
vacuum, it was shown in |^ that the primordial power 
spectrum modified by trans-Planckian effects is given by 



(1) 

where e is the standard inflationary slow-roll param- 
eter related to the inflaton potential ¥{(/)) by e = 
(167r)-iM|(F'/P")2, Mp = is the reduced 

Planck mass, ^ ~ 4 x 10~*y/e/j, 7 = A/Mp parametrizes 
the scale where trans-Planckian effects enter, and Po{k) 
is the primordial power spectrum predicted by standard 
inflationary theory, typically given by the scale-free form 
Po{k) cx fc"=~^, where Ug is the scalar spectral index. 
Finally, fc„ corresponds to the largest scales measurable 
in the CMB. The effect is seen to be a modulation of 
Po(fc)j Sind the natural question to ask is whether this is 
detectable. 

Based on simulated data for the Planck satellite 
it was concluded in Q that the modulation should be 
detectable in the CMB data from this upcoming mission. 
However, this conclusion was based on a Fisher matrix 
analysis, not on a full likelihood analysis. In fact, as 
we will show, the full likelihood analysis reveals that the 
issue is more complicated. The likelihood oscillates in 
parameter space, and one can be misled if one looks only 
at its local behaviour. 



III. OBSERVATIONAL DATA 

A. Cosmic microwave background 

The CMB temperature fluctuations are conveniently 
described in terms of the spherical harmonics power spec- 
trum 
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Since Thomson scattering polarizes light there are ad- 
ditional power spectra coming from the polarization 
anisotropics. The polarization can be divided into a curl- 
free {E) and a curl (_B) component, yielding four indepen- 
dent power spectra: Ct.i, Ge.i, Gb,i and the temperature 
-E-polarization cross-correlation Cte,i- 

The WMAP experiment have reported data on Ct.i 
and Cte,u as described in Ref. 9, 10, 11, 12, IJJ 

We have performed the likelihood analysis using the 
prescription given by the WMAP collaboration which 
includes the correlation between different Q's [gl llollTll 
^1^3]. Foreground contamination has already been sub- 
tracted from their published data. 



In parts of the data analysis we also add other CMB 
data from the compilation by Wang et al. jl^ which 
includes data at high I. Altogether this data set has 28 
data points. 



B. Large scale structure 

The 2dF Galaxy Redshift Survey (2dFGRS) [H has 
measured the redshifts of more than 230 000 galaxies 
with a median redshift of Zm ~ 0.11. An initial estimate 
of the convolved, redshift-space power spectrum of the 
2dFGRS has been determined [ig for a sample of 160 
000 redshifts. On scales 0.02 < k < 0.15/i Mpc"^ the 
data are robust and the shape of the power spectrum is 
not affected by redshift-space or nonlinear effects, though 
the amplitude is increased by redshift-space distortions. 
A potential complication is the fact that the galaxy power 
spectrum may be biased with respect to the matter power 
spectrum, i.e. light does not trace mass exactly at all 
scales. This is often parametrised by introducing a bias 
factor 



h\k) 



Pm(fc)' 
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where Pg{k) is the power spectrum of the galaxies, and 
Pmik) is the matter power spectrum. However, we re- 
strict our analysis of the 2dFGRS power spectrum to 
scales k < 0.15 hMpc^^ where the power spectrum 
is well described by linear theory. On these scales, 
two different analyses have demonstrated that the 2dF- 
GRS power spectrum is consistent with linear, scale- 
independent bias [l7l[l8l |. Thus, the shape of the galaxy 
power spectrum can be used straightforwardly to con- 
strain the shape of the matter power spectrum. However, 
when looking for modulations or other features in the pri- 
mordial power spectrum using the 2dFGRS, one should 
bear in mind that what is measured is the convolution of 
the true galaxy power spectrum with the 2dFGRS win- 
dow function W 0|, 



,(k) « J Pg(k 



q)\Wk{q)\'d'q 



(5) 



and it was found in 19] that this convolution washes out 
any features in the primordial power spectrum for k < 
0.1ft- Mpc^^. However, combining the 2dFGRS power 
spectrum with CMB data breaks parameter degeneracies 
that are present if each dataset is analysed separately, 
and therefore a combination of large-scale structure and 
CMB data gives tighter constraints on the primordial 
power spectrum than the CMB alone. 



C. Likelihood analysis 

For calculating the theoretical CMB and matter power 
spectra we use the publicly available CMBFAST package 
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[20j. As the set of cosmological parameters we choose 
fim, the matter density, flf,, the baryon density, Hq, the 
Hubble parameter, r, the optical depth to reionization, 
Q, the normalization of the CMB power spectrum, b, 
the bias parameter, and finally e and ^. We restrict the 
analysis to geometrically flat models flm + = 1- 

In principle one might include even more parameters 
in the analysis, such as r, the tensor to scalar ratio of 
primordial fluctuations. However, r is most likely so 
close to zero that only future high precision experiments 
may be able to measure it. The same is true for other 
additional parameters. Small deviations from slow-roll 
during inflation can show up as a logarithmic correc- 
tion to the simple power-law spectrum predicted by slow- 
roll. |2U_p2, 23] or additional rclativistic energy density 
in uKm, 2LM,2R.-,MM.M could be present . How- 
ever, there is no evidence of any such effect in the present 
data and therefore we restrict the analysis to the "mini- 
mal" standard cosmological model. 

In this full numerical likelihood analysis we use the 
free parameters discussed above with certain priors deter- 
mined from cosmological observations other than CMB 
and LSS. In flat models the matter density is restricted 
by observations of Type la supernovae to be ^Im = 
0.28 ± 0.14 [ail- The current estimated range for fife/i^ 
from BBN is nth"^ = 0.020 ± 0.002 fsl, and finally the 
HST Hubble key pro ject found the value of Hq to be 
72 ± 8 kms~^ Mpc~^ |33- The marginalisation over pa- 
rameters other than e and ^ was performed using a sim- 
ulated annealing procedure ■ 



IV. RESULTS 

In this section we present the results of the likelihood 
analysis of our selection of current cosmological data sets. 
First of all, in figure ^ we show 68 % and 95 % exclusion 
limits in the ^-e parameter space (all other parameters 
have been marginalised over) for the WMAP TT and TE 
power spectra only. Figure [21 shows how the contraints 
change when the compilation of pre- WMAP CMB data 
from J and in figure 13 the 2dFGRS power spectrum 
data have been included as well. 

From these figures it is clear that there is no preference 
in the data for an oscillating power spectrum, at any 
significant level. The e = ^ = point is always within 
the 1(7 allowed range. 

For the CMB-only data the current excluded range is 
for ^ > 0.1 and 0.1 ^ e ^ 1. ^ is a measure of the am- 
plitude of oscillations, and the bound on this parameter 
comes from the fact that small amplitude oscillations are 
not visible, given the precision of current experiments. 
Since e roughly determines the period of the power spec- 
trum oscillations it is clear why only a very limited range 
is excluded. For small values of e the power spectrum 
does not go through a whole oscillation within the visi- 
ble range of /c-values. On the other hand, for very large 
e the oscillation period becomes small compared to the 
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FIG. 1: 68% and 95% confidence exclusion plot of the param- 
eters e and ^ for WMAP data only. The straight lines are 
isocontours for the parameter 7. 



width of the power spectrum window function. 

In the case where LSS data is included the conclusion 
is very similar, except for the presence of an excluded 
region beyond ^ = 1 at low e. The reason for this sim- 
ply is that when 1^ > 1 the power spectrum can take 
negative values for some k. Since the amplitude fitting 
of the spectrum is done my multiplying with a positive, 
real number this means that no good fit can possibly be 
obtained when starting from an underlying power spec- 
trum with negative values. Of course it also reflects the 
transition to the unphysical region of parameter space 
since the power spectrum should be a positive definite 
quantity. 



V. SIMULATED, HIGH-PRECISION DATA 
SETS 

A. Simulated Planck data 

We have simulated a data set from the future Planck 
mission using the following very simple prescription: We 
assume it to be cosmic variance (as opposed to fore- 
ground) limited up to some maximum Z-value, Zmax, 
which we take to be 2000. For the sake of simplicity 
we shall work only with the temperature power spec- 
trum, Ct,/, in the present analysis. In fact the Planck 
detectors will be able to measure polarization as well as 
temperature anisotropies, but our simplification of using 
only Ct.i will not have a significant qualitative impact 
on our conclusions. For a cosmic variance limited full 
sky experiment the uncertainty in the measurement of a 
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FIG. 2: 68% and 95% confidence exclusion plot of the pa- 
rameters e and ^ for the WMAP data set combined with the 
pre-WMAP data compilation from Wang et al. The straight 
lines are isocontours for the parameter 7. 



rameters can be measured is probably on the optimistic 
side. 



B. Simulated SDSS data 

The Sloan Digital Sky Survey (SDSS) ^ aims at mea- 
suring the redshift of approximately 1 million galaxies, 
and from this the galaxy power spectrum can be obtained 
with unprecedented accuracy. We will consider how this 
improves our ability to constrain trans-Planckian physics, 
assuming that the bias of galaxies with respect to dark 
matter is simple and scale- independent. 

For a fixed survey strategy, the survey volume and the 
number density of galaxies in the redshift sample set a 
lower bound on the uncertainty in the estimated power 
spectrum. In the limit of a perfectly spherical volume- 
limited sample, the uncertainty in the estimated power 
per mode is roughly 



p{k) ~ V 14 



S{k) 



p{k) 



(7) 
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FIG. 3: 68% and 95% confidence exclusion plot of the param- 
eters e and ^ with the previous CMB data plus additional data 
from the 2dF galaxy survey. The straight lines are isocontours 
for the parameter 7. 



given Ci is simply 



Ci 



21 + 1 



(6) 



It should be noted that taking the data to be cosmic 
variance limited corresponds to the best possible case. 
In reality it is likely that foreground effects will be signif- 
icant, especially at high I (see Ref. js^)- Therefore, our 
estimate of the precision with which the oscillation pa- 



[391 where P{k) and P{k) are the true and estimated 
power spectra, S{k) = 1/n is the shot noise power 
for a mean galaxy density n (we will use tT = 2 x 
10^^ h'^ Mpc~^ ), Vc — {2ttY /Vs is the coherence volume 
in the Fourier domain for a survey with volume Vs^ and 
it is assumed that we average the power estimates over a 
shell in Fourier space with volume 14 ~ 47rA:^AA:. That 
is, we average the power over all angles, and over bins 
with AA: > 27r/i?, where R is the survey depth, which we 
wiU take to be 500 /iMpc"^ 

In constructing a mock window function we follow [39| , 
so the window function in Fourier space is the Fourier 
transform of the product of the survey mask with the 
redshift selection function (p(r). Since we are anyway 
interested in constructing an optimistic estimate of the 
window function, we assume that the survey geometry is 
simple and has full coverage of 2tt steradians, 40] and 
adopt a magnitude hmit Bj < 18.9. 

With these optimistic assumptions, the effect of the 
window function is less dramatic than in the case of the 
2dFGRS. In figure 01 we show the ratio of the power spec- 
trum with the standard, scale- invariant primordial Poik) 
to a modulated one with e = 0.03, ^ = 0.023, before 
(solid line) and after (dashed line) convolution with our 
mock window function. As can be seen, the main effect is 
some smoothing of large-scale power, but the modulation 
is still clearly visible. 



C. The CMB window function 

CMB data do not directly measure the underlying pri- 
mordial power spectrum of fiuctuations. Rather, they 
measure the spectrum folded with a transfer function in 



5 



1.04 



1.02 - 



2 1.00 



0.96 
10 





^ — \ 






// 




// 




// 




// 




// 




// 




' / 





0.010 



10 



k (hMpc ) 



FIG. 4: Ratio of modulated to umodulated power spectrum, 
before (solid line) and after (dashed line) convolution with 
the mock SDSS window function. 



the following sense 
Ci 



dk 



(8) 



where A^(fc) is the transfer function taken at the present, 
T = tq, where r is conformal time. Following the line-of- 
sight approach pioneered by Seljak and Zaldariagga 
this transfer function can be written as 



A(fc) 



dTS{k,T)jl[k{T - To)], 



(9) 



where 5* is a source function, calculated from the Boltz- 
mann equation, and ji{x) is a spherical Bessel function. 
However, in order to get a very rough idea about the 
effective window function wi{k) = Af{k)/k of CMB we 
approximate S with a constant to obtain 
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FIG. 5: Ratio of modulated to umodulated power spectrum, 
both the exact result (solid line) and the approximation Eq. 
(16) (dashed line). 



we can write the change in the CMB power spectrum, 
AC I as 



AC, Ii7- ^^dk 



C, 



Po(fc) 



//ro fe3 



dk 



(12) 



wi(k) cx 



fc3 



for fcro < / 
for /ctq I 



(10) 



This simple equation shows several things: (a) A fea- 
ture at some specific wavenumber k = has the great- 
est impact on the CMB spectrum at I* ~ fc*To. For 
a flat, matter dominated universe, tq = Hq/2, yielding 
U ~ 2fc,/_ffo- (b) The CMB windows function is quite 
broad, and narrow features in P{k) are accordingly diffi- 
cult to detect. 

Starting from the power spectrum given in Eq. Q with 
P{k) = Po(fc)(l + AP(fc)) and 



(11) 
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Since kg 1/tq this can be recast in the relatively simple 
form 



AG 1 / , -.-3 • 



r 



dq. (13) 



In the limit where | ^^^2^^^ I 3> 1 one can use the approxi- 
mation 



= — I COS 

Ci e 



21' 



(14) 



In figure we show the quantity AC; / Ci for the same 
two sets of parameters as in Fig. 2 of Ref. Q, i.e. (e = 
0.01,7 = 0.01) and (e = 0.01,7 = 0.003). 
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D. Results 

In figures and [7| we show the result of a hkehhood 
analysis for simulated Planck and SDSS data. The un- 
derlying model used was a standard concordance ACDM 
model {rirn = 0.3, = 0.7, nth^ = 0.024, Us = 1, and 
Ho — 70 kms~^ Mpc~^), with the addition of a modu- 
lated spectrum with parameters, e = ^ = 0.01. From the 
analysis of the simulated Planck data set we find that 
the best fit model for e = ^ = 0.01 has = 2063.1 
for 1/ = 1992 degrees of freedom. This is roughly within 
the expected 1 a interval of = ± V^i' ~ 1992 ± 63. 
However, the actual best fit point of the analysis was at 
e = 0.022,^ = 0.053 which has = 2046.5, seemingly 
ruling out the correct model at more than 99.9% confi- 
dence. 

The problem with this analysis is that both models 
are decent fits within the expectations, but that the ex- 
act value of the likelihood function is extremely sensi- 
tive to e and ^. This is partly due to the fact that the 
phase of the spectrum oscillation at a given fc-value is 
strongly dependent on e and ^. It should also be noted 
that any grid-based likelihood calculation is likely to fail 
when faced with this type of likelihood function, and 
that only stochastic algorithms like Markov Chain Monte 
Carlo (MCMC) or simulated annealing are likely 
to work. 

If, instead of using the standard prescription of assign- 
ing likelihood according to Ax^ compared with the global 
best fit, one uses a robust estimation such as goodness- 
of-fit the plot becomes very different. Figure |S1 shows 
the contours when taking the Icr region to be 1992 ± 63, 
i.e. the expected la interval for the distribution. In 
this case the data is unable to discriminate the under- 
lying model from one without any oscillation, but the 
huge number of non-connected regions in the likelihood 
fit disappear. 

This shows very clearly that the likelihood estimator 
is in principle very powerful, but that it is not robust. 
Using a more robust estimator significantly reduces the 
sensitivity, but on the other hand also eliminates the pos- 
sibility of unphysical biasing of the results. 

From figure |S1 it can also be seen that the exclusion 
region has the same topology as those using the present 
WMAP data (figures 1 and 2), and for the same reasons. 

The analysis is almost equivalent when including the 
mock SDSS data. It should, however, be noted that the 
likelihood contours in figures IHI and [7| are very different, 
showing again the sensitivity of Ax^ to e and ^. Figures 
1^ and [TUI illustrate the degeneracy between models with 
widely different e and ^. 

VI. DISCUSSION 

We have performed a detailed analysis of how a loga- 
rithmic oscillation of the primordial fluctuation spectrum 
influences CMB and LSS observations. The motivation 
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FIG. 6: 68% and 95% confidence exclusion plot of the param- 
eters e and ^ for the simulated Planck data only. The straight 
lines are isocontours for the parameter 7. 
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FIG. 7: 68% and 95% confidence exclusion plot of the param- 
eters e and ^ for the simulated Planck data combined with 
the simulated SDSS data. The straight lines are isocontours 
for the parameter 7. 



for this study was that trans-Planckian effects could lead 
to exactly such a modulation of the spectrum, but in fact 
our analysis of present and simulated future data sets is 
more general. 

The analysis of present data shows that there is no evi- 
dence at a significant level for any oscillating behaviour of 
the power spectrum, but that only a very narrow region 
of parameter space is currently excluded. This region is 
not close to the expected values of the parameters e and 
7 for plausible models of Planck scale physics, such as 
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FIG. 8: 68% and 95% confidence exclusion plot of tiie param- 
eters e and ^ wlien tlie Icr region is taken to be 1992 ± 63, 
tlie expected la interval for the distribution. The straight 
lines are isocontours for the parameter 7. 



the Horava-Witten model. 

In order to study the possibihty of detecting modula- 
tions with future data sets, such as those from the Planck 
satellite or the SDSS survey, wc have constructed mock 
data sets and performed likelihood analyses on them. 
This analysis showed that the likelihood function is ex- 
tremely sensitive to the parameters e and ^, and that 
the likelihood function exhibits a large number of dis- 
tinct local minima. This makes any analysis based on a 
Ax^ approach very difhcult and any result is likely to be 
biased by unphysical effects. On the other hand, more 
robust methods, such as relying on a calculation of 
alone without relying on A^^ as compared to the global 
best fit, are much less sensitive. We have shown for one 
specific case, where parameters were chosen to be close 
to optimistic estimates of the Horava-Witten model, that 
using a standard likelihood analysis yields a biased pa- 
rameter estimate, whereas a robust method is unable to 
distinguish the underlying model from one with no oscil- 
lations. 

Our conclusion is that CMB and LSS data are in prin- 
ciple very sensitive to modulations in the underlying pri- 
mordial power spectrum, but that in practise it is likely 
to be extremely difhcult to make a positive detection of 
the small amplitude oscillations predicted, even with fu- 
ture high precision data. 
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FIG. 9: The CMB power spectra for the best fits to the sim- 
ulated data based on the reference model e = ^ = 0.01 (blue 
line). The green line has e = 10~^'^, ^ = 10"'^ '^, while the red 
line has e = 10"^ \ i = 10"^'^'^. The bottom panel shows the 
residuals with respect to the reference model 
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FIG. 10: The galaxy power spectra for the best fits to the 
simulated data based on the reference model e = ^ = 0.01 
(blue line). The green line has e = lO"^'^, ^ = 10"" ■^ while 
the red line has e = 10"^ S ^ = 10"^'^^ The bottom panel 
shows the residuals with respect to the reference model. 
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